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Abstract 



We study the inclusive two jet triply differential cross section d^a/dETdr]idr]2 at 
Fermilab energies. Different r]i and r]2 pseudorapidity regions are directly related to 
both the parton level matrix elements and the parton densities at leading order. We 
present the next-to-leading order [©(af )] corrections and show that the shape of the dis- 
tribution at fixed transverse energy Et is a particularly powerful tool for constraining 
the parton distributions at small to moderate x values. We investigate the renormal- 
isation/factorisation scale uncertainty present in the normalisation and shape of the 
distribution at next-to-leading order. We discuss specific slices of the distribution, the 
same-side/opposite side ratio and the signed pseudorapidity distribution, in detail and 
compare them with preliminary experimental data. 



1 Introduction 



Dijet production in hadron collisions occurs when two partons from the incident hadrons 
undergo a hard pointlike interaction and scatter at relatively large angles. The two-jet cross 
section depends on both the non-perturbative probability of finding a particular parton inside 
the parent hadron and the dynamics of the hard scattering. By examining kinematic regions 
where the parton densities are well known, we can probe the pointlike strong-interaction 
cross section. One example of this is the shape of the angular distribution of the jets in their 
centre of mass frame. Recent data provide clear evidence that a running coupling constant, 
as given by QCD, is needed to describe the data []T| and that the next-to-leading order QCD 
predictions are in good agreement with the data. An alternative approach is to use the 
theoretical description of the hard scattering to extract the distribution of partons in the 
proton from the data. This is particularly interesting since gluon scattering plays a very 
important role in two jet production, and it may be possible to probe the gluon density in 
a more direct way than is possible in deeply inelastic scattering or in Drell-Yan processes. 

The inclusive two-jet cross section can be described in terms of variables most suited 
to the geometry of the detector; the transverse energy of the leading jet, Et = Eti, and 
the pseudorapidities of the two leading jets, rji and 772. Recently, the DO collaboration has 
presented a preliminary measurement of d^a/dETdr)id'r]2 as a function of t]i and 772 at fixed 
Et- This seemingly complicated three-dimensional quantity contains all the information 
available from two jet events. In particular, at leading order, t]i and r]2 are directly related 
to the parton momentum fractions Xi, X2, 



so that a measurement of the triply differential cross section d^a / dETdr]idri2 at fixed E^ 
corresponds to a measurement of d'^a / dxidx2- Although the overall normalisation of cross 
sections is uncertain in perturbative QCD, one might hope that the shape of this distribution 
is well-predicted and that it can be used to discriminate between different parton densities. 
Inclusion of the next-to-leading order corrections enhances the reliability of the calculation 
for both shape and normalisation. 

Beyond leading order, however, Eq. || is no longer satisfied and the parton momentum 
fractions are only approximately determined by the transverse energies and pseudorapidities 
of the jets. A three-dimensional plot may obscure some of the desired physics, so both 
the CDF 1^ and DO [|l], Q collaborations have focused on particular slices of the general 
distribution. The CDF collaboration has examined the ratio of cross sections for same-side 
events (jji ~ 772) to opposite-side events (?7i ~ —'U2) for different Et bins. This reduces the 
normalization uncertainty and enhances the small x region, x ~ AE^ / s. At small transverse 




(exp(r/i) +exp(r72)) , 
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(exp(-?7i) + exp(-r72)) , 
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energies, this distribution can reliably discriminate between singular {xg{x) ~ at small 
x) and non-singular xg{x) ~ x° behavior of the gluon distribution |^. The DO collaboration 
has taken slices in rji which contain information over the whole x range, AE^/s < x < 1, 
but are more sensitive to the overall normalisation. In this paper, we study the shape 
of the distribution over the whole 'r]i-r]2 plane. First, we examine the triply-differential 
distribution at lowest order. We discuss how the available phase space grows as the number 
of final state partons increases and relate different rji and ri2 regions to both the parton 
level matrix elements and the parton densities (section 2). The full next-to-leading order 
triply-differential cross section is presented in section 3. We show that it is sensitive to 
the parton density functions and indicate how the shape depends on the renormalisation 
(and factorisation) scale. Our results are applied to the CDF same-side/opposite-side and 
DO signed distributions in sections 4 and 5. Finally, our main findings are summarised in 
section 6. 



2 The 0{a^^) triply differential two jet cross section 

The lowest order cross section is given by, 

dETdr]idr]2 8^ Ej. cosh 77* 

where fi{x,fiF) ii = g,<l,<l) represents the density of parton i in the proton at factorisation 
scale /ii? and |A^jjP is the lowest order squared matrix element for ij 2 partons summed 
and averaged over initial and final state spins and colours. The strong couping constant as 
is evaluated at the renormalisation scale, /ir. The parton level cross section is insensitive 
to Lorentz boosts and |A1jjP therefore depends only on the parton pseudorapidity in the 
parton-parton rest frame, rj* = (r^i— ?72)/2. To understand how the cross section is distributed 
over the rji — 772 plane, we recall the 'single effective subprocess approximation' 0. In 
this approximation all parton-parton scattering cross sections are taken to be equal but 
are weighted by colour arguments. Thus the gluon-gluon, quark-gluon and quark-quark 

subprocesses are in the ratio 1 : | : (|) so that. 



dETdr]idri2 Stt ' ' Ef. cosh 77 

where F{x,fi) is the 'single effective parton density', 

4 

F{x, fi) = g{x, /i) + - X! (^(^' /^) + ^(^' /^)) • (4) 

y q 

A rough indication of how the physical cross section depends on r]i and rj2 can be obtained 
by studying the parton-parton luminosity, XiF(xi,/i) X2F{x2, and the squared matrix 
elements, \M.gg\'^/ cosh"^?]*, for gluon-gluon scattering separately. 
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Figure 1: The phase space boundary in the rji — 772 plane at leading-order (solid) and next- 
to-leading order (dotted) for Et — 50 GeV and = 1800 GeV. The dashed line separates 
the 'small' x - 'large' x and 'large' x - 'large' x regions. In region I, either xi or X2 is less 
than Xt, while in region II, both xi and X2 are bigger than xt- Region III is only permitted 
at next-to-leading order. 



However, to orient ourselves in the rji — rj2 plane, we first focus on the allowed phase 
space in terms of rji and 772. At lowest order, the jet pseudorapidities are directly related to 
the parton fractions via Eq. 1. Since the momentum fraction cannot exceed unity, we find. 




Xt J \ Xt 

and, 

|,,|<cosh-'(-L), (6) 

where xt = IEtI \fs and x\ < X1X2 < 1. This boundary is shown in Fig. 1 for Et = 50 GeV 
and y/s — 1800 GeV. For the opposite-side cross section, rji ~ —772, the parton fractions are 
roughly equal so that in the top left and bottom right corners of the allowed phase space, 
Xi ^ X2 ^ 1. On the other hand, in the bottom left and top right corners, corresponding to 
same-side events with rji ri2, the parton fractions are maximally different, Xi x^, X2 — > 1 
and vice versa. The dashed boundary separating region I and II makes a nominal division 
of the phase space according to whether both parton fractions are 'large' or one parton 
fraction is 'small'. In region II, xi and X2 > xt while in region I, either xi < xt or X2 < xt- 
The corresponding axes for the pseudorapidities of the two jet system in the laboratory, 
^boost = (^1 + ^2)72, and of the jet in the jet-jet center of mass frame, rj* = {rji — 772)72, are 
related by a rotation of 45°. 
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Figure 2: Contours of (a) constant r]i and (b) constant rj* in the X1—X2 plane for Et = 50 GeV 
and y/s = 1800 GeV. 

At lowest order in perturbative QCD, each point in the r/i — 772 plane is uniquely related 
to the parton momentum fractions Xi and X2 via Eq. 1. To give an idea of how different 
X values are spread over the allowed pseudorapidities, we show contours of fixed rji in the 
Xi — X2 plane in Fig. 2. At fixed rji, the smallest xi value occurs when ri2 is a minimum, 
while T]2 is a maximum for the smallest X2 value. An alternative way of looking at the xi —X2 
plane using the rj* and r^boost variables is shown in Fig. 2b. For a given r]*, varying r^boost 
over its allowed range takes {xi,X2) = (1, xf. cosh^(r7*)) to (xi,X2) = (xf, cosh^(?7*), 1) while 
preserving X1X2 = x'^ cosh"^ {r]*) . These contours are particularly useful since the parton 
matrix elements depend only on rj* and not on T^boost- 

Although the charged parton distributions have been probed directly over a wide range 
of parton momentum fractions x and scales in deeply inelastic scattering, the gluon 
density is rather poorly known. Direct photon data from WA70 determine the shape of 
the gluon in the x ~ 0.3 — 0.4 region, however, the gluon density at other x values is only 
constrained by the momentum sum rule. To explore the sensitivity of the triply differential 
cross section, we choose parton density functions with contrasting small-x behaviours; the 
improved MRSD„ and MRSDq distributions of ref. ^ for which xg{x) behave as x~^'^ and 
x° respectively at small x and Q^. The low- a; behaviour of measured at HERA is better 
fitted by an x~^'^ growth as parameterised by the MRSA distributions ||^, however the 
range of predictions from the MRSD_ and MRSDq distributions indicate where the triply 
differential cross section is sensitive to the small-x parton distributions. It is worth noting 
that because of the momentum sum rule, a parton density that is relatively large at x ~ 10~^ 
must be relatively small at x ~ few x 10~^. This is demonstrated in Fig. 3a where we show 
the ratio of the 'single effective parton density' for the MRSDq and MRSD_ parton density 
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Figure 3: (a) The ratio of the 'single effective parton density' of Eq. 4 for the MRSDq and 
MRSD_ distributions compared to the MRSA parameterisation at // = 50 GeV. (b) The 
ratio of the gluon and quark parton densities in the MRSDq distribution compared to the 
MRSA parameterisation at the same scale. 

functions relative to that of the MRSA parameterisation. The hierachy evident at small x 
is reversed in the moderate x range, where the MRSDq density is 15% larger than that for 
the MRSD_ density functions. 

At X ~ 1, the MRSD distributions are both larger than MRSA. This is primarily be- 
cause the up and down valence distributions arc fitted separately in the more recent MRS 
parameterisations rather than the up and the (up plus down) valence combination. In any 
event, the distribution of partons inside the proton with x > 0.5 is very small and poorly 
constrained by data. 

If we examine the gluon and the quark densities separately (the latter summed over all 
quark and anti-quark flavors), shown in Fig. 3b, then we find that other than at small x, 
the quark densities are quite well determined by present-day data: different sets are quite 
similar. In contrast, it is the gluon densities that are poorly determined: different sets are 
substantially different even at intermediate x. As suggested by the 'single effective parton 
density,' and as we shall see in greater detail, the gluon densities are sufficiently important 
to jet production in hadron-hadron scattering to cause substantial variations in predictions 
dependent on the densities at moderate x. 

As can be seen from Eq. 2, the cross section is proportional to the product of structure 
functions. To get a feehng for how this product varies. Fig. 4 shows the parton-parton 
luminosity in the single effective subprocess approximation as a function of 77boost for different 
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Figure 4: The parton-parton luminosity for the MRSDq and MRSD_ parton densities in the 
'single effective subprocess approximation' as a function of 77boost for I??*! — 0, 1 and 2 and 
/X = 50 GeV. 



|77*| values. This corresponds to diagonal strips across the r)i — 772 plane. As expected, the 
largest luminosity occurs when Xi and X2 are equally small, 77* ~ 7/boost = 0. Once again, the 
MRSDq luminosity is approximately 20% larger than that for MRSD_ . As either |?7boost| 
or increases, the luminosity decreases rapidly. However the falloff is more rapid with 
increasing \r]*\ than with increasing Ir^boostl- 

As mentioned earlier, the parton-parton subprocess scattering matrix elements are inde- 
pendent of r/boost- One consequence is that for fixed 77*, the only variation of the cross section 
comes from the variation of the parton densities as jyboost runs over the allowed kinematic 
range. The lowest order squared matrix elements for gg — > gg scattering are given by, 

\M,f ^ Qtt^ (4cosh^(r7*)-l)=^ 
cosh^(r/*) 8 cosh^(r;*) ' ^ ^ 

These matrix elements are plotted in Fig. 5 as a function of f]*. At small \ri*\ the matrix 
elements grow rapidly until |?7*| ~ 2 where the matrix elements saturate. This behaviour 
complements the parton-parton luminosity which is largest at rj* — 0. 

By multiplying the parton-parton luminosity with the squared matrix elements (along 
with the overall factor a'j.{n)/8nEj. ) we obtain the physical cross section. Because different 
parton distributions dominate for different momentum fractions (and hence rji and 772 values), 
we expect the shape of the triply-differential cross section to be sensitive to the parton 
densities. This is illustrated in Fig. 6, where we show the leading order prediction for the 
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Figure 5: The lowest order squared matrix elements cosh^(?7*) for gg gg scattering 

as a function of rj*. 

ratio, 

^ _ Ie^- dEr (^.^^(MRSDo) - ;,^^(MRSD_)) 

/l^^-- dEr , (MRSD,) ' ^ ^ 

J-C-Tmin ^ (IEt dr]i dr]2 ^ ' 

in the transverse energy range 45 GeV < E^ < 55 GeV evaluated at /i = Et^. In addition to 
the sharp cutoff marking the boundary of the allowed phase space, the excess of MRSDq over 
MRSD_ at small r/i and ri2 and the depletion at large |?7i| ~ \ri2\ are seen clearly. Over the 
whole rji — 772 plane, the relative cross sections vary by +23% to —5%, with the most sizeable 
effects at r^i ~ 772 ~ where the cross section is largest. At lowest order, Eti = Et2 = Et, 
so that the distribution is symmetric under rji ^ ri2- 



3 The 0{al) triply differential two jet cross section 

At next-to-leading order, some contributions admit three partons into the final state; for 
these, the parton fractions are given by, 

xi,2 = exp{±rii) + -— exp(±r72) + ^ exp(±r73) , (9) 

where Exi and rji {i = 1, ... ,3) describe the transverse energies and pseudorapidities of the 
three partons ordered in decreasing Et- Since the transverse energies of the partons are no 

^The raw cross sections are histogrammed in 0.5 x 0.5 bins in 771 and 772- 
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Figure 6: The leading order prediction for the ratio of the triply differential cross section 
R{LO) defined in Eq. 8 for 45 GeV < Et < 55 GeV and /x = Et- 

longer forced to be equal, |?72| may increase to compensate for having a smaller transverse 
energy, Et2/Eti < 1. The maximum possible values of |?72| occur when Et2 = E^s, 



a j<*<'°s( i j. (10) 

where a = (2 — xti exp(r7i))/xTi and a = (2 — xti exp(— ?7i))/xti. The enlargement of 
phase space is shown in Fig. 1 (region III). We see that the maximum allowed value of 
|?7i| is unchanged at next-to- leading order so that the physical cross section will exhibit a 
rather sharp cutoff as \rii\ increases. On the other hand there will be a more gradual fall 
off in the cross section as |?72| increases. Indeed, adding more partons into the final state 
further increases the allowed r/2 range corresponding to the production of more and more 
soft partons. 

To compute the next-to-leading order cross section, we use an 0{a^) Monte Carlo pro- 
gram for one, two and three jet production based on the one-loop 2^2 and the tree 
level 2 — > 3 parton scattering amphtudes [|10|, |ll| described in ref. 0]. This program uses 



the techniques of refs. O, M to cancel the infrared and ultraviolet singularities thereby 



rendering the 2 — 2 and 2 — >^ 3 parton processes finite and amenable to numerical com- 
putation. The parton four momenta are then passed through a jet algorithm to determine 
the one, two and three jet cross sections according to the experimental cuts. Different cuts 
and/or jet algorithms can easily be applied to the parton four-momenta and, in principle, 
any infrared-safe distribution can be computed at 0{a^) 



In order to compare the theory with experiment, we use the parton level equivalent of the 
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Figure 7: The next-to-leading order triple differential distribution for 45 GeV < Et < 
55 GeV and /i = Eti for (a) MRSDq and (b) MRSD_ parton densities. 



standard 'Snowmass' cone algorithm |T5[ with AR = 0.7 and require at least two jets in the 
event. Furthermore, we note that the assignment of which jet is hardest is not infrared safe, 
so that we must symmetrize over the hardest and next-hardest jets (in transverse energy). 
The distributions we shall examine require that the hardest jet lie in a certain 'trigger' range; 
each event will be counted twice if the next-hardest jet also lies in this Et range. The three- 
dimensional cross section evaluated at yU = Eti, where Eti is the transverse energy of the 
hardest jet in the event, is shown in Fig. 7 for 45 GeV < Et < 55 GeV. The extension of 
the phase space to smaller 772 is seen clearly, along with the rather sharp cutoff at \rii\ ~ 3.5. 

Although the cross sections for the two parton densities appear similar, the difference 
between the predictions observed at lowest order is preserved. This is illustrated in Fig. 8, 
where we show the next-to-leading order prediction for the fractional difference ratio defined 
in Eq. 8. As at lowest order, the difference varies between +23% at ?7i ~ 772 ~ and —10% 
at ?7i ~ r72 ~ —2.5. We note that the ratio is most negative when xi ~ Xt and X2 ~ 1. This 
is the region where the singular behaviour of the MRSD_ parton densities dominate over 
the less singular MRSDq distributions 0. 

One indicator of the reliability of perturbation theory is the ratio of next-to-leading order 
to leading order cross sections^ This is shown in Fig. 9 for two slices of the triply differential 
distribution. First we consider the slice rji = and let 772 vary. At large \ri2\, the phase space 
extends beyond the strict leading-order kinematic limit of t]2 < cosh~^(l/xT) = 3.68. As a 

^The normalisation is extremely sensitive to the renormalization scale choice in the leading order cross 
section. 
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Figure 8: The next-to-leading order prediction for the ratio of the triply differential cross 
section R{NLO) defined in Eq. 8 for 45 GeV < Et < 55 GeV and = ^ti- 

consequence, the corrections are large. However, this does not signal the emergence of large 
logarithms which might spoil the applicability of perturbation theory. Rather, the large 
corrections are due to the more restricted phase space available at leading order, that is the 
absence of region III of Fig. 1. Towards the edges of available phase space, the leading-order 
cross section is thus forced artificially to zero, and the ratio goes to infinity. At yet-higher 
order, however, one expects further corrections to be reasonable (until one approaches the 
outer edges of region III or its higher-order analogs). Second, we keep ri2 = fixed and allow 
rji to vary. We see that the corrections for central |?7i| are small, however as the magnitude 
of rji approaches the edge of phase space, the next-to-leading order corrections significantly 
reduce the cross section. In this limit, the second jet is forced to have Et2 ~ Eti and 
the available phase space for soft gluon emission is curtailed. As a consequence, radiative 
corrections lower the cross section close to the edge of phase space. These corrections are 
a result of the appearance of large logarithms, and one does expect perturbation theory to 
behave badly near this edge of phase space. In summary, we see that the triply differential 
cross section is reliably predicted over the whole range of the 771—772 plane with the exception 
of the very large |?7i| slices. 

We have seen how the shape of the two-dimensional distribution is sensitive to the parton 
density functions. However, there is also a dependence on the renormalisation and factori- 
sation scales fiR and //^Qthat could in principle obscure the differences due to the structure 
functions. To get a feeling of how severely the scale uncertainty affects the shape. Fig. 10 

•^Throughout we choose /ii?, = /iF = however, other choices are possible. 
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Figure 9: The ratio of next-to-leading order (NLO) to leading order (LO) predictions as a 
function of 772 for rji = (solid) and of rji for = (dashed) for 45 GeV < Ex < 55 GeV, 
fi = Et and the MRSD_ structure functions. 



shows the next-to-leading order predictions for the ratio, 

for the MRSD_ parton densities and A = 0.5 and 2. Because the absolute magnitude of 
the cross section (which depends on q;s(/x)) is poorly predicted, the prediction for /i = XEti 
has been normalised to the cross section for = E^i at r/i ~ ?72 ~ 0. For A = 0.5, 
c = 0.93, while for A = 2, c = 1.08. We have restricted the pseudorapidity range in the 
plot to Ir^il, |?72| < 2.5 since for higher pseudorapidities the next-to-leading order effects are 
large as discussed above. As a result, for such pseudorapidities, there is a sizeable scale 
variation. However, for central pseudorapidities, the shape is changed by less than 5% which 
is significantly less than the difference between the two representative parton distributions. 

Once the experimental data are available, it should therefore be possible to extract in- 
formation on the density of partons in the proton. In addition to the uncertainty in the 
normalisation of the theoretical predictions, there is a significant uncertainty in the experi- 
mental normalisation as well, due to uncertainties in the luminosity measurement, jet energy 
calibration, jet trigger efficiency, and other aspect^. It therefore makes sense to allow the 
overall normalisation of the theoretical prediction, a^^, to float, so that by varying c the 

"^Part of the experimental uncertainty can be eliminated by normalising with respect to the W cross 
section. 
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Figure 10: The next-to- leading order prediction for the ratio of the triply differential cross 
section R{NLO) as defined in Eq. 11 with (a) A = 0.5 and (b) A = 2 for 45 GeV < Et < 
55 GeV and yU = Eti- 

for 

JETmin ^ \dEj'driidri2 dETdriidri2 

summed over the different r/i, 772 cells is minimised for a given parameterisation of the parton 
densities. Finally, the input parameterisations can be adjusted so that the is further 
reduced. This can be done simultaneously for different slices in transverse energy. An even 
more interesting possibility would be to map out the evolution of the parton densities directly 
by following trajectories of constant (xi,X2) in the rji — 772 plane as a function of Et- 




4 The DO signed distribution 

Recently, the DO collaboration has presented preliminary data Q for a particular slicing of 
the triply differential distribution — the so-called signed pseudorapidity distribution. This 
amounts to taking two strips of the rji — 772 plane for a fixed transverse energy interval and 
combining them in reverse directions. The pseudorapidity of the leading jet is constrained 
to lie in the range |?7i|min < |^i| < I ^1 1 max and the distribution is plotted as a function of 
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Figure 11: The signed pseudorapidity distribution for (a) 45 GeV < Et < 55 GeV, 0.0 

< Ir^il < 0.5 and (b) 55 GeV < Et < 65 GeV, 2.0 < < 2.5 at both LO and NLO. The 
prehminary experimental results from are also shown. The factorisation scale is chosen 
to be the transverse energy of the hardest jet, jU = Eti- 

1 772 1 sign (771772), 

da _ 1 /"'^^"''"'^^ [^/'''^'"''"'(i / '''^'""'"d '^^^ \ 

(i 1 7/2 1 sign (771 772) AETJEr^in ^ 2Ar]i \J\,jiUn ^ dETdr]idr]2 ^-ImUax ^ dETdr]idr]2 ) '' 

(13) 

where sign(?7i772) = —1 if ?7i and 772 have opposite sign and +1 if they have the same sign. 
Positive values of I772I sign (771772) correspond to same-side dijet events, while negative values 
are associated with opposite-side events. In principle, both strips contain equal information, 
but combining them serves to reduce the statistical error. Once again, we must sum over 
the hardest and next-hardest jet in order to ensure that this distribution is infrared-safe. 

In the currently available data, DO has examined two slices in transverse energy, 45 GeV 

< < 55 GeV and 55 GeV < < 65 GeV, and two strips in 7/1, 0.0 < |7/i| < 0.5 and 
2.0 < |7/i| < 2.5. As more data become available from the current Tevatron run, this analysis 
can be extended to cover a larger range of Et and rj. 

We first fix 7/1 to lie in the central pseudorapidity slice, 0.0 < |7/i| < 0.5 and examine the 
pseudorapidity of the second jet for the transverse energy interval 45 GeV < Et < 55 GeV. 
This strip includes the 771 ~ 7/2 ~ region that is sensitive to parton densities at a; ~ 0.05. 
We also consider the slice at larger 771, 2.0 < \rii\ < 2.5 but for a slightly higher transverse 
energy interval, 55 GeV < Et < 65 GeV. The predictions for these distributions for both 
MRSDq and MRSD_ parton densities are shown in Fig. 11 with the preliminary data from the 
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Figure 12: The (a) parton-parton luminosity for the MRSD_ parton densities and (b) squared 
matrix elements in the 'single effective subprocess approximation' for t]i = 0.5 and Et = 
50 GeV (solid) and rji = 2.5 and Et = 60 GeV (dashed). The factorisation scale is chosen 
to be 50 GeV. 

DO collaboration [jl|0. We see a clear asymmetry favouring smaller values of 1 772! sign (771772 )• 
In other words, for rji ~ 0, the opposite-side cross section (negative \r]2\sign{riiri2)) peaks 
away from 772 ~ 0, while the same-side cross section (positive |?72|sign(?7i?72)) monotonically 
decreases. A similar effect has been observed in the same-side/opposite-side cross section 
measured by CDF 0. 

This is is due to an interplay between the parton-parton luminosity and the matrix 
elements. Fig. 12 shows the parton-parton luminosity in the 'single effective subprocess 
approximation' and the gg gg matrix elements of Eq. 7 as a function of 1 172 1 sign (171772). 
We see that the maximum of the parton-parton luminosity occurs at 1 772 1 sign (7/1 772) ~ 0.25 
and 0.75 for 771 = 0.5 and 2.5 respectively, while the minimum of the matrix elements always 
lies at |772|sign(77i772) = |77i|. The net effect of the shift in the peak of the parton-parton 
luminosity to positive 1 772 1 sign (771 772) combined with the shift of the minimum of the matrix 
elements to larger values of |772|sign(77i772) is an enhancement of the cross section at negative 
1 772 1 sign (771 772) and a depletion at positive 1 7/2 1 sign (771772), clearly visible as an asymmetry in 
Fig. 12. 

As suggested by the solid line in Fig. 9, the next-to-leading order corrections reduce the 
cross section uniformly by about 10% until the kinematic limit on 772 from the lowest order 

^We have divided the data by a factor of two to account for the size of the pseudorapidity interval. 
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Figure 13: The ratio of ncxt-to- leading order predictions for the signed distribution for the 
MRSD_ parton densities evaluated at = XEti for A = 0.5 (solid), A = 2 (dashed) and 
A = 4 (dotdashed) relative to that for fi = Eti ^ot (a) 45 GcV < Et < 55 GeV and 
0.0 < Ir^il < 0.5 and (b) 55 GeV < Et < Q5 GeV and 2.0 < l^^il < 2.5. The ratio of next- 
to-lcading order predictions for the MRSDq and MRSD_ parton densities with /j, = Eti is 
shown dotted line. 

process is approached. We also see that the difference between the MRSDq and MRSD_ 
predictions is about 23% at I772I ~ as expected from Fig. 8. 

At larger pseudorapidities, the next-to-leading order predictions give a much better de- 
scription of the data than at leading order. The preliminary data appear to favour the 
MRSDq parameterisation at the x values probed here, x ~ 0.05. However, the errors are 
still large, and as mentioned at the end of section 3, there are significant uncertainties in the 
overall normalisation of the experimental data. 

As discussed in the previous sections, there is also an uncertainty in the normalisation 
of the theoretical cross section due to the choice of renormalisation and factorisation scales. 
This is particularly evident for the signed distribution since the lowest order cross section 
is proportional to o^^iir). Even at next-to- leading order, the overall normalisation is still 
uncertain. However, one would expect that the shape of the distribution is relatively insen- 
sitive to varying hr- This is illustrated in Fig. 13, which shows the ratio of next-to-leading 
order predictions for different scales relative to the next-to-leading order MRSD_ prediction 
for 11 = Eti- 

As expected, the normalisation of the small \r)2\ region is quite sensitive to the choice of 
scale, however, the shape of the distribution for small \rj2\ is essentially unchanged for the 
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central rji slice, 0.0 < \rii\ < 0.5. For large pseudorapidities, \ri2\ > 3, there is of course 
a large scale dependence. This is a consequence of exceeding the lowest-order kinematic 
limit on m; in this region, an 0{al) calculation such as the one performed here is in fact a 
leading-order one. As a contrast, the ratio of the next-to-leading order predictions for the 
MRSDq and MRSD_ parton densities at /x = Eti are also shown. In addition to a sizeable 
change in the normalisation, the shape of the distribution around |?72| ~ is also changed. 
It remains an experimental question as to whether this difference in shape can be detected. 



5 The CDF same-side over opposite-side ratio 

The interpretation of DO signed distribution measurement hinges strongly on the absolute 
normalisation of the cross section. Thus one would need to know the jet energy correction 
well before one can constrain the parton density functions. To circumvent this problem the 
CDF collaboration has considered a ratio, that of same-side (SS) to opposite-side (OS) cross 
sections 0. For the same-side cross section, both jets have roughly the same pseudorapidity, 
while in the opposite-side cross section the jets are required to have roughly equal, but 
opposite pseudorapidities. One then forms this ratio, as a function of the pseudorapidity in 
several tranverse-energy slices. In a realistic experimental analysis, the pseudorapidities and 
transverse energies will be binned so that. 



F = / / , (14) 

. rq+Ari /---q+Ari rETmi^^ d^fj 

(Tosim^ = / dT]i dT]2 dEr . (15) 

^ETmin<ET<ETma^ Jr^-Ari J-Ti-Ari JErmin dETdriidri2 



From these cross sections we form the SS/OS ratio. 



R. 



SS/OS 



Ermin < Et < Eti, 



ETmin<ET<ETn 



(^osiv) 



(16) 



ETmin<ET<ETn 



with the advantage that a large part of the experimental and theoretical uncertainties cancel. 
However, most of the dependence on the parton densities in the central region where r]i ~ 
~ is also removed. As we saw in the previous sections, this is exactly the region where 
there can be a strong dependence on the parton density functions. Nevertheless, we can still 
study the behaviour of the gluon density at small x by examining the SS/OS ratio at large 
pseudorapidity. Since the x values probed are much smaller than in the signed pseudorapidity 
distribution, typically x ~ AE^/s rather than x ~ lE^j \fs., studying this ratio is to a large 
extent complementary to studying the DO signed pseudorapidity distribution. 

In the preliminary CDF measurement 0], the jet transverse energy was chosen to lie in 
four separate bins, 27 GeV < < 60 GeV, 60 GeV < < 80 GeV, 80 GeV < Et < 
110 GeV and 110 GeV < Et < 350 GeV and the pseudorapidity interval to be Ar/ = 0.2. 
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Figure 14: The leading order parton fractions probed by the SS/OS cross section ratio for 
Et = 15 GeV. 

It may prove possible to extend the analysis to smaller transverse energy and include a fifth 
bin, 15 GeV < Et < 27 GeV. As for the triply-differential and signed rapidity distributions, 
we must symmetrize explicitly over the leading and ncxt-to-leading transverse energy jets 
in order to ensure that the sided cross sections are infrared safe. In addition, in order to 
suppress events with three or more hard jets, an azimuthal angle cut between the two leading 
jets of TT - 0.7 < A0 < TT + 0.7 is applied. 

To get a feeling for the range of parton fractions probed by this particular cross section. 
Fig. 14 shows the three different momentum fractions (using the leading order definition of 
Eq. 1) as a function of the pseudorapidity for the smallest accessible jet transverse energy 
Et = 15 GeV. We see that it is possible to probe parton fractions for x values as small as 
3 X 10"^ <x<l. 

The next-to-leading order predictions for the SS/OS ratio for the very low transverse 
energy bin, 15 GeV < Et < 27 GeV, are shown in Fig. 15. Because very low x values 
are encountered, the different parton density functions give a broad range of predictions. At 
present, no data are available for this particular transverse energy range. However, even with 
relatively large experimental uncertainties one can still easily discriminate between different 
parton densities. This makes the very small Et bin particularly interesting. 

The region around rj = 3 has the largest sensitivity to different parton density functions, 
corresponding to a smallest parton fraction of 8 x 10~^. At larger pseudorapidities (and 
therefore smaller x values) there is a severe phase space supression and the cross section 
(and event rate) decreases rapidly. 
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Figure 15: The next-to-leading order predictions for the SS/OS ratio evaluated at = Eti 
for the MRSD_ (sohd), MRSDq (dotted) and MRSA (dashed) parton distributions as a 
function of rj for the smallest Et bin accessible to the CDF collaboration, 15 GeV < Et < 
27 GeV. 

In the other four transverse energy bins, the CDF collaboration has published prelimi- 
nary data which we display along with next-to-leading order predictions for three different 
parton density sets in Fig. 16. Because the experimental results are preliminary, we should 
be careful in drawing conclusions from these results for the time being. These results do 
however prefigure discriminatory powers which should emerge as the experimental errors 
shrink with the inclusion of more data from the current Tevatron run. Even with the cur- 
rent uncertainties, the smallest E^ bin suggests a deficiency in the density of MRSDq partons 
round x ~ 10~^, and in contrast to the signed rapidity distribution discussed in section 4, 
favors the MRSD_ set of structure functions. 

For the higher transverse energy bins higher x values arc sampled and there is not much 
difference between the different parton density functions. As an example, the highest trans- 
verse energy bin is mostly sensitive to parton momentum fractions around 0.1, where the 
parton density functions are, in principle, tightly constrained. Of course, one should still 
compare the data with theory. That different parton densities give the same results docs 
not guarantee they are correct; it merely indicates that they either use the same data as a 
constraint or use the same assumptions to derive the individual parton densities. 

For this distribution as well, we should also consider the uncertainty in the theoreti- 
cal predictions arising from renormalisation and factorisation scale dependence present in 
perturbative QCD calculations. To study this. Fig. 17 shows the SS/OS ratio for the 27 
GeV < £'t < 60 GeV transverse energy bin at pseudorapidity 77 = 2.6 as a function of 



15GeV<ET<27 GeV 



■ NLO, MRSD- 
NLO, MRSDO 
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Figure 16: The next-to- leading order (NLO) predictions for the SS/OS ratio evaluated at 
/i = Eti for the MRSD_ (solid), MRSDq (dotted) and MRSA (dashed) parton distributions 
as a function of t] for the four transverse energy bins, 27 GeV < Et < 60 GeV, 60 GeV 
< Et < 80 GeV, 80 GeV < < 110 GeV and 110 GeV < Et < 350 GeV. Also shown is 
the preliminary CDF data of ref. 



= I^R = I^F- This phase space point offers the highest discriminatory power with the 
currently available data. As an illustration we have chosen the MRSD_ set which seems 
to be favoured by these data. The reference scale is the transverse energy of the highest 
Et jet in the event, /i = Eti, as in the previous sections. The leading-order prediction 
does not depend on the strong coupling constant. As the SS and OS cross sections probe 
the parton densities at different momentum fractions x, it does depend on the factorization 
scale, but this gives rise only to a trivial scale dependence: the SS to OS ratio rises lin- 
early with log{fi/ETi). At next-to-leading order, the prediction does depend on the couphng 
constant, and the overall dependence is less trivial. It so happens that the reference scale, 
fi = Eti, coincides with the minimum in the variation of the ratio with respect to scale. 
Either increasing or decreasing the scale results in an increase of the SS/OS ratio. Varying 
the scale by a factor of two around the reference scale, we get a feeling for the theoretical 
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H/Et 

Figure 17: The renormalisation/factorisation scale dependence of the SS/OS ratio for rj = 2.6 
in the 27 GeV < Et < 60 GeV bin using the MRSD_ structure functions. The data point 
is taken from ref. 0. 

uncertainty. The largest variation comes from reducing /i and increases the cross section 
by approximately 20%. Due to this uncertainty one cannot really discriminate between the 
MRSD_ distributions with xg{x) ~ and the more recent MRSA fits with xg{x) ~ x~^'^. 
However, it should be possible to exclude the MRSDq distributions once the final CDF data 
are published. To discriminate between MRSD_ and MRSA, it will be necessary to make a 
measurement in the 15 GeV < Et < 27 GeV transverse energy bin. 

Fig. 18 shows the ratio of next-to-leading to leading order predictions (the i^-factor) for 
27 GeV < Et <Q0 GeV and 80 GeV < Et < 110 GeV using the MRSD_ parton densities 
and fi = Eti- The shape of the SS/OS ratio is basically unchanged between —2 < rj < 2. 
For larger pseudorapidities we get a rapid change in the ii'-factor. This is again mainly due 
to the fact that the leading order cross section is quickly forced to zero by the kinematic 
constraints on 2 — 2 scattering, and does not indicate the presence of large logarithms which 
might spoil the applicability of perturbation theory. We have parameterised the ii'-factor as 
an even polynomial in rj, 

K{r])=A + B7f + Cri^ + Dri^ + Erf. (17) 

The fitted constants A . . . E for all five transverse energy intervals are given in Table 1. We 
see that for this double ratio, the corrections are extremely small in the central region. 
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Figure 18: The ratio of next-to-leading to leading order predictions for the SS/OS ratio with 
// = Eti and MRSD_ parton densities in the five transverse energy intervals. 

6 Conclusions 

In this paper we have made a detailed study of the two-jet cross section in hadron-hadron 
collisions for a given range of jet transverse energy as a function of the pseudorapidities of 
the two jets. This distribution is particularly sensitive to the parton density functions at 
small {x ~ few x 10^^) to intermediate {x ~ few x 10^^) parton momentum fractions for 
jet transverse energies accessible to the CDF and DO experiments at Fermilab. For example, 
parton distributions that are relatively large at small x (and therefore constrained to be rel- 
atively small at larger x by the momentum sum rule) lead to a relative enhancement of the 
same-side two-jet cross section at large pseudorapidities {r]i ~ 772 ^ 0) and a relative deple- 



Pseudorapidity bin 


ABODE 


15 GeV <Et<21 GeV 
27 GcV < < 60 GeV 
60 GcV < Et < 80 GcV 
80 GcV < Et < 110 GcV 
110 GeV <Et< 350 GeV 


1.0197 0.0062 0.0226 -0.0025 8.166e-5 

0.9992 -0.019 0.0453 -0.0094 5.592e-4 
0.9653 0.072 0.0326 -0.0208 0.0027 
0.9869 0.0891 -0.0332 0.0057 6.427e-4 
1.0121 -0.1974 0.4985 -0.2743 0.0457 



Table 1: The parameterisation of the i^T-factor for MRSD_ in the 5 transverse energy bins. 
The fitted formula is K{r]) — A-\- Brf -\- Cr]'^ -\- Drf -\- Erf and is intended to be used for 
returning the X-f actor for the bin, given the center of the bin (that is 77 = 0.2, 0.6, 1.0, ...) 
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tion for central jets, r^i ~ ?]2 ~ 0. We have studied the next-to-leading QCD corrections to 
this distribution in some detaiL To discuss the properties of the cross section it is convenient 
to divide the rji — rj2 plane into two regions, the central region where the pseudorapidity of 
the two jets is less than the lowest order kinematic limit, \r]i\, \r]2\ < cosh~"'^(l/a;x) and the 
foward region where the rapidity of one of the jets approaches or exceeds the lowest order 
boundary. In the central region, 

i) the next-to-leading order corrections are small and perturbation theory works well. 

ii) the scale uncertainty in the overall normalisation is reduced. Varying by a factor of 2 

about fjL — Eti changes the next-to-leading order prediction by 0(8%). 

iii) the scale uncertainty in the shape of the distribution is quite small. For the same 
variation of ji, the relative bin-to-bin correction is less than C(2%) — see Fig. 10. 

On the other hand, in the forward region, 

i) the next-to-lcading order corrections are important and improve the agreement with 

experimental data. 

ii) there is a considerable scale uncertainty because the corrections calculated in a 'next- 

to-leading' order program corrections are effectively lowest order. 

Because the central pseudorapidity region is both sensitive to the parton densities and stable 
to higher order corrections, the triply differential distribution offers an excellent chance to 
gain extra information about the distribution of partons in the proton. The current and 
future runs of the Tevatron at Fermilab should yield copious quantities of two-jet events and 
once precise experimental data are available, it should be possible to make a determination 
of the gluon density at small and intermediate x values. However, in addition to the uncer- 
tainty in the normalisation of the theoretical predictions, there is also an uncertainty in the 
experimental normalisation due to uncertainties in the luminosity measurement, in the jet 
energy calibration, and other effects. The information on the parton densities lies more in 
the shape than the overall normalisation, and one way of determining the parton densities 
is to allow the overall normalisation of the theoretical prediction, cr^^, to float, so that by 
varying c the for 



summed over the different 771, r]2 cells in the central pseudorapidity region is minimised for 

a given paramctcrisation of the parton densities. By adjusting the input paramctcrisations, 
the be further reduced as the predicted shape becomes closer and closer to that 

observed in the data. Of course, data can be taken for many slices in transverse energy. 
This allows the possibility of following the evolution of the parton densities directly by 
following trajectories of constant (xi, x^) in the 771 — 772 plane as a function of Et- 
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At present only preliminary data are available for particular slices of the r]i — r]2 plane. 
The signed distribution (sect. 4) presented by the DO collaboration is primarily sensitive to 
intermediate x values. As we may expect from the preceding discussion, the shape is well 
predicted, but the normalisation is uncertain. Nevertheless, the preliminary data appear 
to favour the MRSDq paramcterisations (Fig. 11). On the other hand, the same-side to 
opposite-side cross section ratio presented by the CDF collaboration probes much smaller x 
values. Once again, the shape is relatively unchanged by including the next-to-leading order 
corrections. As shown in Fig. 16, the preliminary data (again with large errors) appear to 
favour the more singular MRSD_ paramcterisation. Of course, with the current experimental 
data sample, no definitive conclusion can be drawn. However, if these tentative observations 
are accurate, it would imply that the data favour the parton density with the largest density 
of gluons at both x ~ few x 10~^ and x ~ few x 10~^. In other words, there are more 
gluons present in the small and intermediate x regions than expected from the momentum 
sum rule, which may in turn suggest that the density of gluons is not as well determined 
by direct photon data (WA70) as previously thought. Data from the current Tevatron run 
should help to provide an answer to this puzzle. 
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